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The development of a prophylactic vaccine that targets human papillomaviruses (HFV) 6, 
11 , 16, and 18 to prevent cervical cancer has increased interest in the ethnic and 
geographical distributions of HPV genotypes. We investigated HFV prevalence and type 
distribution by restriction fragment mass polymorphism (RFMP) testing a total of 60,775 
specimens (aged 18-79 yr, median 44) taken from liquid-based cytology. Overall HFV 
positive rate of total patients was 34.2%. Among the positive patients, 87.7% was single 
type infections, and 12.3% was multiple HFV types. HFV-16 was the most prevalent 
genotype observed in 2,307 (26.0%), followed by type 52 in 2,269 (25.5%), type 58 in 
1 ,090 (12.3%), type 18 in 633 (7.1%), type 56 in 436 (4.9%). The pattern of high risk- 
HFV positive rate according to age showed U-shape with a peak in HFV prevalence among 
women less than 30 yr of age, and a second peak among the older females aged 70 to 79 
yr. The leading four high-risk HFV genotypes were HFV-16, HFV-52, HFV-58, and HFV-18 
in descending order. In conclusion, this study provides the most representative prevalence 
and type-specific distribution of HFV among Korean women, and demonstrates that the 
epidemiology of HFV infection is different from that of other regions of the world. 
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INTRODUCTION 

Infection with human papillomavirus (HPV) is considered to 
be a pre-requisite for the development of cervical cancer and 
is associated with neoplastic precursor lesions (1). In Korean 
women, the cervical cancer is the fifth frequent cancer after 
breast cancer, stomach cancer, colon cancer and thyroid can- 
cer. To date more than 100 HPV genotypes have been identified 
and at least 50 are known to infect the female anogenital tract 
(2). Among these thirteen mucosotropic HPVs (types 16, 18, 31, 
33, 35, 39, 45, 5, 52, 56, 58, 59, and 68) have been classified as 
class I high risk carcinogens to human beings (3). Because the 
distribution and prevalence of HPV vary by geographic region 
and the immunity conferred by vaccines is type-specific, the 
need for HPV genotyping in routine screening population is in- 
creasingly recognized. For example, although HPV-16 is the 
most prevalent type worldwide, some genotypes such as type 



52 and 58 are rare in Western countries, whereas they are rela- 
tively prevalent in Asian populations (4-6). Hence, an accurate 
assessment of the regional, community-based distribution of 
HPV genotypes is extremely important for the prevention of 
cervical cancer and for public hygiene management. The US 
FDA-approved Hybrid Capture II HPV DNA test (HC II) (Digene 
Corporation, Gaithersburg, MD, USA), which can detect 13 car- 
cinogenic HPV types, is the HPV DNA detection method most 
commonly used. Unfortunately, this cocktail detection method 
does not identify specific HPV genotypes. The restriction frag- 
mentation mass polymorphism (RFMP) assay (GeneMetrix Co., 
Yongin, Korea), which utilizes PCR to amplify the HPVL1 gene 
followed by enzyme restriction and mass measurement using 
MALDI-TOF, offers a method to detect 38 types of HPV. This 
new method has manifested sensitivity and reliability compa- 
rable to HC II (7-9). 

The goal of this study was to establish and investigate an age- 
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specific HPV prevalence, genotype distribution and extent of 
multiple infections in Korean women to provide the most rep- 
resentative epidemiologic data for the clinical prediction of the 
outcome of HPV infection, future prevention strategies, and the 
development of the efficient multivalent HPV vaccine. 

MATERIALS AND METHODS 
Study population 

The residual samples after liquid-based cytology (LBC) tests 
from 60,775 patients aged 18-79 yr (median age 44-yr-olds) from 
September 2006 to September 2011 were used for HPV geno- 
types using RFMP assay. After cytology, residual cells in the LBC 
sample were centrifuged at 3,500 rpm for 10 min and stored as 
split cellular pellets of 200 uL at -70°C prior to nucleic acid ex- 
traction. 

DNA extraction 

QIAampl DNA Blood Mini Kits (Qiagen, Valencia, CA, USA) 
were used to extract cervical cell DNA according to the manufac- 
turer's protocol. Nucleic acid was stored at -70°C prior to HPV 
amplification and genotyping. 

HPV DNA amplification and genotyping 

The presence and genotype of HPV in all samples were tested 
using RFMP method. The details of the RFMP assay protocol 
were described previously (8). Briefly, 4 aL of DNA were ampli- 
fied with PGMY09/1 1 primers, comprising of two non degener- 
ate pools of LI consensus primers. Second round primer pairs 
comprised a sense primer RFMP specific to bases 6584 to 6603 
(5 '-GCMCAGGGHCAYAAGGATGAA TGG-3 ') and an antisense 
primer RFMP specific to bases 6657 to 6626 (5'-GTACTDCK- 
DGTRGTATCHACMACG GATGTAACAAA-3 ')• The 5-nucleo- 
tide sequence (GGATG) embedded in the primers introduced a 
Fokl site (a neoschizomer of BtsCTj in the amplicon. Restriction 
enzyme digestion of PCR products was performed by mixing the 
PCR reaction mixtures with 10 uL of buffer containing 50 mM 
potassium acetate, 20 mM Tris-acetate, 10 mM magnesium ac- 
etate, 1 mM dithiothreitol and 1 unit of Fokl and BtsCl. The re- 



action mixtures were incubated at 37°C for 1 hr. The resulting 
digest was purified by vacuum filtration through 96-well Oasis® 
aElution Plates (Waters, Tokyo, Japan). The desalted reaction 
mixtures were resuspended with matrix solution containing 
15 mg/mL 3-hydroxypicolinic acid, 0.023 M ammonium citrate 
and 12% acetonitrile, and spotted in 3-uL volumes on a polished 
MTP AnchorChip™ plate (BrukerDaltonics, Bremen, Germa- 
ny). Mass spectra were acquired with the aid of installed soft- 
ware (flexcontrol 3.0) on a Microflex linear MALDI-TOF mass 
spectrometer (BrukerDaltonics, Bremen, Germany).This meth- 
od was able to detect 38 HPV genotypes (6, 11, 16, 18, 26,30,31, 
33, 34, 35, 39, 40, 42, 43, 44, 45, 51, 52, 53, 54, 55, 56, 58, 59, 61, 
62, 66, 68, 70, 72, 73, 74, 81, 82, 83, 84, 89, and 90). 

Statistical analysis 

The data were analyzed by an SPSS 12.0 statistical package for 
Windows. Pearson's chi-squared test was used to evaluate the 
significant differences between the designated groups. All tests 
were two-sided, and a P value < 0.05 was considered statistical- 
ly significant. 

Ethics statement 

This study was approved by the institutional review board of 
Green Cross Reference Laboratory of Korea (GCRL-IRB-2012001). 
Informed consent was exempted by the board due to the retro- 
spective design and anonymous data. 



Table 1. Overall HPV prevalence stratified by age 

HPV genotyping result by RFMP assay 



Age (yr) 


No. positive (%) 


No. negative (%) 


Total No. 


18-29 


2,143(49.9) 


2,154 (50.1) 


4,297 


30-39 


6,093 (36.0) 


10,814(64.0) 


16,907 


40-49 


6,352 (30.7) 


14,309 (69.3) 


20,661 


50-59 


4,365 (32.2) 


9,187 (67.8) 


13,552 


60-69 


1,195 (33.1) 


2,417(66.9) 


3,612 


70-79 


639 (36.6) 


1,107 (63.4) 


1,746 


Total 


20,787 (34.2) 


39,988 (65.8) 


60,775(100) 



RFMP, restriction fragmentation mass polymorphism; HPV, human papillomavirus. 



Table 2. LR-HPV and HR-HPV positive rate stratified by age 







HPV (+) 








Age (yr) 


LR-HPV infection 


HR-HPV infection 




HPV (-) 


Total 




(single + multiple) 


(single + multiple) 


Total 






18-29 


972 


1,171 


2,143 


2,154 


4,297 


30-39 


2,698 


3,395 


6,093 


10,814 


16,907 


40-49 


3,126 


3,226 


6,352 


14,309 


20,661 


50-59 


2,443 


1,922 


4,365 


9,187 


13,552 


60-69 


654 


541 


1,195 


2,417 


3,612 


70-79 


266 


373 


639 


1,107 


1,746 


Total 


10,159 (16.7%) 


10,628 (17.5%) 


20,787 (34.2%) 


39,988 (65.8%) 


60,775 (100%) 



LR-HPV, low risk-HPV; HR-HPV, high risk-HPV. 
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Table 3. Single and multiple type infections stratified by age 



HPV (+) 







Cjnnlp twnp 
on iyic LyptJ 






Mi iltinlp tvnp 




HPV(+) Total 




LR-HPV (%) 


HR-HPV (%) 


Total 


LR-HPV (%) 


HR-HPV (%) 


Total 


18-29 


881 (49.7) 


893 (50.3) 


1,774 


91 (24.7) 


278 (75.3) 


369 


2,143 


30-39 


2,489 (46.6) 


2,856 (53.4) 


5,345 


209 (27.9) 


539(72.1) 


748 


6,093 


40-49 


2,91 1 (51 .3) 


2,766 (48.7) 


5,677 


215(31.9) 


460 (68.1) 


675 


6,352 


50-59 


2,249 (58.5) 


1,597 (41.5) 


3,846 


194(37.4) 


325 (62.6) 


519 


4,365 


60-69 


587 (56.3) 


456 (43.7) 


1,043 


67 (44.1) 


85 (55.9) 


152 


1,195 


70-79 


232 (42.3) 


317(57.7) 


549 


34 (37.8) 


56 (62.2) 


90 


639 


Total 


9,349 (51 .3) 


8,885 (48.7) 


18,234 (87.7) 


810(31.7) 


1 ,743 (68.3) 


2,553 (12.3) 


20,787 (100%) 



LR-HPV, low risk-HPV; HR-HPV, high risk-HPV. 
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Overall HPV positive rate (%) Table 4. Number of genotypes in multiple infections stratified by age 



Multiple infections of HPV 
positive patients (%) 



-— o- 



■o 



18-29 30-39 
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Age groups (yr) 
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Fig. 1. Overall HPV positive rate and multiple infection rate of HPV positive patients 
stratified by age. 

RESULTS 



HPV prevalence 

Of the 60,775 Korean women aged between 18 and 79 yr, HPV 
DNA was detected in 20,787 (34.2%); 10,628 (17.5%) patients 
were infected with high-risk HPV genotypes and 10,159 (16.7%) 
patients with low-risk HPV genotypes (Table 1, 2). Among HPV 
positive patients of 20,878, 18,234 (87.7%) patients had single 
infection and 2,553 (12.3%) had multiple infection (Table 3). Age- 
specific HPV infection prevalence peaked at 49.9% in 18- to 29-yr- 
olds, and declined to 30.7% in 40- to 49-yr-olds, after which HPV 
prevalence steadily increased to 36.6% in 70- to 79-yr-olds (Table 
1, Fig. 1). In the multiple infections of 2,553, infections with two 
different genotypes were 2,347 (91.9%) and with three genotypes 
206 (8.1%) (Table 4). The number of genotypes was not signifi- 
candy different among the different age groups (P = 0.26, Table 
4). Age-specific prevalence curve of multiple infections in HPV- 
positive patients revealed a very similar pattern to that of the 
overall HPV positivity stratified by age (Fig. 1). 

Age-specific HPV genotype distribution 

We considered HPV genotypes 16, 18, 31, 33, 35, 39, 45, 51, 52, 
56, 58, 59, and 68 as primary carcinogenic (high-risk) types and 





Number of genotypes (%)* 




Age (yr) 


2 


3 


Total 


18-29 


335 (90.8) 


34 (9.2) 


369 (14.4) 


30-39 


700 (93.7) 


47 (6.3) 


747 (29.3) 


40-49 


620 (91 .7) 


56 (8.3) 


676 (26.5) 


50-59 


469 (90.4) 


50 (9.6) 


519(20.3) 


60-69 


142(93.4) 


10(6.6) 


152(6.0) 


70-79 


81 (90.0) 


9(10.0) 


90 (3.5) 


Total 


2,347(91.9) 


206 (8.1) 


2,553 (100.0) 


*X, P=0.26. 








Table 5. The distribution of high-risk HPV genotypes of HPV positive patients 


HR-HPV 


Number of patients (%) 


Total (%) 


Single type (%) 


Multiple type (%) 


16 


2,307 (26.0) 


411 (23.6) 


2,718(25.6) 


18 


633 (7.1) 


160(9.2) 


793 (7.5) 


31 


383 (4.3) 


144(8.3) 


527 (5.0) 


33 


208 (2.3) 


67 (3.8) 


275 (2.6) 


35 


312(3.5) 


60 (3.4) 


372 (3.5) 


39 


125 (1.4) 


29 (1 .7) 


154 (1.4) 


45 


203 (2.3) 


54(3.1) 


257 (2.4) 


51 


405 (4.6) 


102(5.9) 


507 (4.8) 


52 


2,269 (25.5) 


407 (23.4) 


2,676 (25.2) 


56 


436 (4.9) 


112(6.4) 


548 (5.2) 


58 


1,090 (12.3) 


128(7.3) 


1,218(11.5) 


59 


113(1.3) 


20 (1.1) 


133 (1.3) 


68 


401 (4.5) 


49 (2.8) 


450 (4.2) 


Total 


8,885 (100.0) 


1,743 (100.0) 


10,628 (100.0) 



HR-HPV, high-risk human papillomavirus. 

all others as non-oncogenic (low-risk) types (3). 

Among the high-risk HPV (HR-HPV) infections, HPV 16 (25.6 
%) and HPV 52 (25.2%) were the most commonly detected types, 
followed by HPV 58, 18, 56, 31, 51, 68, 35, 33, 45, 39, and 59 at 
11.5%, 7.5%, 5.2%, 5.0%, 4.8%, 4.2%, 3.5%, 2.6%, 2.4%, 1.4%, and 
1.3%, respectively (Table 5). The leading four HR-HPV genotypes 
stratified by age in descending order were HPV type 52, 16, 58, 
and 18 in 18- to 29-yr-olds, and 40- to 59-yr-olds, HPV 16, 52, 
58, and 18 in 30- to 39-yr-olds, and 60- to 69-yr-olds, and HPV 
16, 52, 58, 56, and 68 in 70- to 79-yr-olds (Fig. 2). 
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HR-HPV genotype distribution in aged 18 to 29 yr 
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HR-HPV genotype distribution in aged 30 to 39 yr 
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HR-HPV genotype distribution in aged 40 to 49 yr 
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HR-HPV genotype distribution in aged 50 to 59 yr 
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4.2 


11.8 








^^^m 5.5 






3.4 












0.2 






3.8 






M 2.3 






3.6 






^^^^^m 8.0 






i 



0.0 



10.0 20.0 
% 



30.0 



68 


5.4 








59 
58 


0.0 


1 14.2 






56 
52 


5.4 




23.3 




51 


3.8 








45 


■ 1.6 








39 


1 0.3 








35 


H 2.5 








33 


3.5 








31 


M 2.2 








18 


M 2.8 








16 


III 


0.0 10.0 


20.0 


30.0 


40.0 



Fig. 2. The distribution of high-risk HPV genotypes according to age groups. 



DISCUSSION 

HPV is one of the most common sexually transmitted infections 



worldwide (10). The estimated risk of HPV is notably high over 
a woman's lifetime (over 80%), however, most women who ac- 
quire HPV infection do not develop more serious high-grade 
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cervical neoplasia or invasive cancer because most infections 
are transient (11). 

Our cross-sectional study was performed on patients referred 
to the gynecology office seeking for an opportunistic cancer 
screening with liquid based cytology. Population-based data 
for HPV-type distribution is a prerequisite to the development 
of new HPV vaccine and cancer screening tests, and to assess 
the effect and benefits of vaccination. The prevalence of cervi- 
cal infection with human papillomavirus (HPV) in women var- 
ies geographically (12). This study represents, to our knowledge, 
the largest study of HPV prevalence, including HPV-DNA prev- 
alence and a type specific distribution stratified by age in 60,775 
women in Korea. Therefore, the present study provides a more 
closely represented HPV prevalence of the general population 
in Korean women for vaccine development and future prophy- 
lactic programs in Korea. The cytology results of all samples were 
composed of 84% of normal cytology, and 16% of abnormal cy- 
tology (data not shown). 

Aside from geographical differences, variation may also be 
the result of detection methods employed and varied specimen 
type and amount. Therefore, precise and accurate HPV DNA 
test depends on the selection of a proper HPV genotyping meth- 
od. Although the FDA-approved Hybrid Capture II (HC II) HPV 
DNA test is the most widely used HPV DNA detection method, 
it does not have the capability to identify a specific HPV geno- 
type. Furthermore, previous studies have shown that the HC II 
HPV test cross-reacts with at least 15 HPV genotypes not includ- 
ed in its current high-risk probe cocktail set (13, 14). In this study, 
HPV genotypes were identified by using RFMP assay which is 
comparable to HC II with regard to detection of high-risk HPV 
infection, but is also able to identify 38 HPV genotypes in a sin- 
gle assay, which offers high analytical and clinical sensitivity and 
the advantage of reliable detection of multiple HPV infections, 
and it showed a good correlation with direct sequencing data 
(7, 9). Several studies on the prevalence of overall HPV infection 
in Korean women have reported rates ranging from 7.2% to 44.8%, 
and the prevalence of HR-HPV infection ranging from 7.2% to 
34.3% (15-20). Because the prevalence of HPV infection varies 
by country, region within country, detection method, sample 
size and specimen type, the prevalence of HPV in Korean women 
shows in wide range according to various articles. As this study 
was conducted with a very large number of samples of 60,775, 
using validated precise and accurate detection method, the prev- 
alence of HPV infection observed on this investigation might be 
the most accurate. Overall, HPV infection rate of our study was 
34.2%, and high-risk HPV infection positivity was 17.5%. Those 
results are lower than the previous study in which overall HPV 
prevalence accounted for 44.8% and high-risk HPV infection 
rate was 34.3% (19). Age-specific HPV prevalence varied con- 
siderably across geographical regions. Across all geographical 
regions, observed HPV prevalence was strongly associated with 



age, and the curves of HPV prevalence according to age were 
characterized by a U-shape with a relatively higher HPV preva- 
lence in younger and older ages in Africa, Central/South Amer- 
ica, Europe, and North America (21). This study also revealed 
such U-shape curve, which peaked at women between 18- and 
29-yr-olds, decline women between the age of 40- and 49-yr- 
olds, and then steadily increased to women between the age of 
70- and 79-yr-olds (Fig. 1). The U-shaped curves of HPV by age 
may potentially be explained by newly acquired HPV infections 
and reactivation of latent HPV infections in the older women 
(22, 23). 

Multiple infections among HPV positive cases were accounted 
for 12.3%, which was lower than previous study which revealed 
20.9% (14). Several studies reported that the risk for developing 
the tissue abnormalities or lesions that typically precede cervi- 
cal cancer is much higher for women infected with multiple gen- 
otypes of the human papilllomavirus. Therefore, incorporating 
diagnoses of multiple HPV infections into the prediction out- 
comes of HPV infections is considered important (24, 25). The 
prevalence of multiple HPV infections in HPV positive cases has 
shown a worldwide geographical variation ranging from about 
9% to 50% in the European countries (26, 27). This variation is 
also due to different sample size, population, detection meth- 
ods and the type of the samples. 

Distribution of HPV genotypes varies across geographic ar- 
eas (12, 27, 28). Several studies demonstrated that the most com- 
mon type is HPV 16 followed by 18 in Europe, Central and South 
America, by HPV 52 and HPV-58 in Asia, by HPV-53 and HPV 
52 in North America (6, 10, 27, 29). These differences in HPV 
type distribution in countries and regions may be related to dif- 
ferent sexual habits and migrations of people (27, 30). Therefore, 
accurate information on the distribution of HPV genotypes based 
on the population is very important for both primary cervical 
cancer screening and prophylactic vaccination policy decision 
making. 

Data on the distribution of HPV genotypes in Korean women 
remain controversial. Cho et al. showed that the three major 
HPV types were HPV 16, 18, and 33 (17). Shin et al. showed that 
the most common HPV types in southern Korean women were 
70, 16, and 33 (20). Hwang et al. reported that three major HPV 
types were HPV- 16, HPV-52, and HPV-58 in descending order, 
which is the same result of our study. The study by Hwang et al. 
is based on relatively large sample size of 2,470 (19). 

Our study showed that the five most prevalent high-risk HPV 
types were Type 16, 52, 58, 18, and 56 in descending orders. Those 
five HPV types comprised about 76.7%, and four most prevalent 
types comprised 72.0% of all the 13 high risk HPV type infections. 
On the other hand, HPV-16 and 18, which was the most com- 
mon HPV type in the western countries, comprised 34.0% of all 
the 13 high risk HPV type infections. The four most prevalent 
high-risk HPV types were not different across six age groups, 
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except the fourth prevalent high-risk HPV type was HPV-56 and 
68 in women between aged 70- and 79-yr-olds. Therefore, high 
risk HPV types other than HPV-16 and 18 such as genotypes 52, 
58, and 56 must be targeted in the design of diagnostic tests and 
prophylactic vaccines for Korean women. 

HPV genotype prevalence and distribution were different 
compared with the other regions, showing higher frequencies 
of HPV 52, and HPV 58. The rate of multiple genotype infections 
was lower than other countries as 12.3%. The most common 
four HR-HPV genotypes were HPV16, 52, 58, and 18 in descend- 
ing order. There is a peak in HPV prevalence among women less 
than 30 yr of age, with a second peak among older females aged 
70 to 79 yr in Korea. 

In conclusion, this study provides the most representative 
prevalence and type-specific distribution of HPV among Kore- 
an women, and demonstrates that the epidemiology of HPV 
infection is different from that of other world regions. 
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